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Abstract 

We develop numerical methods for approximating Ricci flat metrics on Calabi- 
Yau hypersurfaces in projective spaces. Our approach is based on finding bal- 
anced metrics, and builds on recent theoretical work by Donaldson. We illustrate 
our methods in detail for a one parameter family of quintics. We also suggest 
several ways to extend our results. 
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1 Introduction 

Calabi-Yau manifolds were proven to admit Ricci flat metrics in [1] . Explicit expressions for 
such metrics would have many applications in mathematics and in string compactification; 
however it is widely thought that for compact Calabi-Yau manifolds no closed form expression 
exists, except in trivial cases (tori and orbifolds). Thus it is useful to develop methods for 
constructing and working with approximate Ricci-flat metrics; and for the solutions of related 
equations, such as hermitian Yang-Mills, and the other equations of string compactification. 

One such approximating method is finding the balanced metric on an embedding of the 
manifold into complex projective space, given by the sections of a holomorphic line bundle. 
This method was suggested by Yau in the early 90's, and was proven to work in a fundamental 
paper by Tian [2] . Since then many people have worked on the problem of balanced metrics 
in various contexts [3, 4, 5, 6, 7], just to name a few. A recent result of Donaldson [8] 
shows that this approximation scheme is both mathematically elegant and relatively easy 
to implement numerically. It can also be directly generalized to find approximate hermitian 
Yang-Mills connections on vector bundles [9], which can in turn be used to compute metrics 
on moduli spaces, and the kinetic terms in N = 1 string compactifications [10]. 

In this work we develop numerical methods for approximating Ricci flat metrics on Calabi- 
Yau hypersurfaces based on these ideas. This also supplies a detailed analysis of the numer- 
ical methods used in [9] . We study the effectiveness of our approach in the example of a one 
parameter family of quintics in CP 4 . As we review in section 2, we work with a space of 
approximating metrics parameterized by an N x N hermitian matrix; the balanced metric 
is then the fixed point of the so-called "T map" on this space, defined in Eq. (7). 
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The main computational problem in implementing the T map numerically is the evalua- 
tion of a large number of integrals on the manifold. More precisely, given a Calabi-Yau n-fold 
X, with its corresponding holomorphic n-form Q G Q n,0 (X) = A n (T*X) 1,0 , and volume form 
d/j,n = f2 A f2, one needs to compute integrals of the type 

/ dfm, (1) 



x 

where / : X — > C is a smooth complex valued (but not holomorphic) function. Consequently, 
the heart of the paper (section 3) will be devoted to developing a numerical approximation 
scheme to efficiently and accurately compute these integrals. 

A second technical point, which is very valuable in simplifying these computations, is to 
take advantage of the discrete symmetries of the manifold. We discuss this in section 3.4.3. 

Our explicit numerical results appear in section 4, where we also provide a general dis- 
cussion of the efficiency and accuracy of the algorithm, comparisons with alternatives, and 
suggestions for future work. 

Before we begin, let us briefly set out the problem. Denote the Ricci flat metric on X 
(which is unique given a complex structure and Kahler class) as qrf- We want to propose 
a set of approximating metrics gh parameterized by parameters h, and give a numerical 
procedure to find the "best" approximation to g^p within this set. 

The criteria that a best approximation should satisfy include 

1. Accuracy: we want to minimize the error e = d(gh, 9rf)i where d is some measure of 
the distance between the approximate and true metrics. A simple and natural choice 
for e in the present context is to consider the function 

det u h 

on X, where uo^ is the Kahler form for g^. For a Ricci flat metric, this will be the 
constant function. We then take 1 

min xeX 7] h (x) 

e = l — . (3) 

max x£X Vh{x) 

Of course, one could use other norms, such as \ \r)h ^ J rjh\\ p , or curvature integrals. 

2. Control: we want an explicit bound on the error, 

e {9h,9RF) < (4) 

depending on the parameters of the problem. 

3. Systematic improvement: we would like to have a control parameter k, such that 
by increasing k, we can bring the error estimate e max down to any desired accuracy. 



1 We consider only compact varieties, hence both the minimum and the maximum are attained. 



2 



4. Mathematical naturalness. Our experience with string theory (and more generally 
in mathematics and physics) has been that in exploratory work such as this, rather 
than trying to incorporate all known aspects of a problem and find a "best" solution, 
we can learn far more by studying a well chosen simplification in depth. This favors a 
scheme in which one makes the smallest possible number of arbitrary or ad hoc choices 
not inherent in the original statement of the problem. 

Of course the approximation should be efficiently computable as well. We will comment on 
these various aspects as they arise. 



2 T-map 

In this section we review the construction of balanced metrics, which in turn lead to a 
convenient approximation of the Ricci flat metric on Calabi-Yau threefolds. Our presentation 
is based on [9], which in turn was deeply inspired by [8]. We refer to these papers, and the 
references therein, for more details. 

We start with a holomorphic line bundle £ on a Kahler manifold X, with N global 
sections. This fact is usually phrased as H°(X, C) = C N . Let {s a }^ =l be a basis of this 
vector space, and consider the map 

i : X — P^" 1 i(Z , ...,Z n ) = ( Sl (Z), s 2 (Z), s N (Z)). 

The geometric picture is that each point in our original manifold X (with coordinates Zj) 
corresponds to a point in parameterized by the sections s a . Since choosing a different 
frame for L would produce an overall rescaling s a — » Xs a , the overall scale is undetermined. 
Granting that s±(Z), si{Z), . . . , sn(Z) do not vanish simultaneously, this gives us a map into 

pJV-l 

We want this map to be an embedding, i.e. that distinct points on X map to distinct 
points on P^ -1 , and that tangent vectors are also separated. In general, we can appeal to 
the Kodaira embedding theorem, which asserts that for C ample this will be true for all 
powers C k , starting with some k^. As an example, for non-singular quintics in P 4 , Ox{k) is 
both ample and very ample for all k > 1. 

Next we consider the iV 2 -parameter family of Kahler potentials on P^ -1 , 

K h = log |x> Q ^j = lo sll s ll^ ( 5 ) 

parametrized by the N by N hermitian matrix h. These give rise to a family of Kahler 
metrics on X by restriction, and we will seek a "best" approximation to the Ricci flat metric 
within this space of approximating metrics. Note that it is the inverse of h, h a/3 , that appears 
in (5). The reason for this will become clear shortly. 
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Mathematically, the simplest interpretation of (5) is that it defines a hermitian metric 
on the line bundle £. This is a sesquilinear map from £ <g> £ to smooth functions C°°(X), 
here defined by 



f 8 ,s') = e~ Kh -s-s' ~ 



Notice that a change of frame, which acts on our sections as s a — > Xs a , cancels out in this 
expression. 

This metric allows us to define an inner product between the global sections: 

(sp\s a ) = / ^| dvo\ x . (6) 

Jx \\ s \\h 

Note that this inner product depends on h in a nonlinear way, since h appears in the de- 
nominator, and \\s\\h involves the inverse of h. Here dvolx is a volume form on X, which 
has to be chosen. 

One choice which makes sense for any X is dvolx — detu;^, the standard volume form 
y/det gh- For X a Calabi-Yau variety, we can instead use dvolx = v = Q A Q. The latter is 
significantly simpler, mainly because it is h- independent. Our numerical procedure will use 
a discrete approximation to this volume form. 

The explicit form of the expression in (6) suggests that one studies the map 



T{h) a p = — — — / dvo\ x , (7) 



N f s a s p 
^-^ol(X)jxu^ h 

dubbed the "T-map" , which acts on the space of hermitian matrices in a non-linear way. 

A fixed point of this map, T(h) = h, the pair (h, s a ) is called a balanced embedding, and 
the metric on X associated to the corresponding Kahler potential (5) is called the balanced 
metric. 2 As it turns out, the balanced metric is unique, up to U(N) transformations of 
the basis of sections and rescaling, provided that the manifold in question has no continuous 
symmetries. This is certainly the case for a Calabi-Yau variety, and in particular our quintics. 

It turns out that the T-map is contracting, so the simplest way to find a fixed point of 
the T-map is to iterate it. We have the following 

Theorem 2.1 (see, e.g. ,[13, 6]; and [3, 8] for the u-balanced case) Suppose that Aut(A, L) 
is discrete. If a balanced embedding exists, then, for any initial hermitian matrix Go, as 
r — > oo the sequence T r (G ) converges to a fixed point. 

The importance of balanced metrics stems from a theorem that goes back to Tian [2] 
and Zelditch [5]. Let us consider the sequence of balanced metrics associated to the bundles 



2 In a more precise nomenclature, which we will not use here, the term "balanced embedding" is reserved 
for the definition with dvolx = detoJ/i, while the alternate definition with dvolx = v is referred to as the 
'V-balanced embedding." Also, there are several other equivalent ways of defining the notion of balanced 
metric, which are closely linked to the notion of stability in Mumford's sense, but this utilitarian definition 
will suffice for our purposes. For reviews of this topic we recommend [11, 12]. 
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C k (denned in terms of the Fubini-Study metric oj^ s from Eq. (5)): 

u k = \il{u F k s ). (8) 

The rescaling is made so that the cohomology class of the Kahler form [uik] = C\{C) G 
H 2 (A, Z) is independent of k. With this definition one has 

Theorem 2.2 Suppose Aut(A, £) is discrete and (X,£ k ) is balanced for sufficiently large 
k. If the metrics u>k converge in the C°° norm to some limit as k — > oo, then is a 
Kahler metric in the class c\(C) with constant scalar curvature. 

The constant value of the scalar curvature is determined by c\{X). In particular, for 
c\{X) = the scalar curvature is zero. Thus, the balanced metrics Uk, in the large k limit, 
converge to the Ricci flat metric. Furthermore, in this case the 1 jk convergence is enhanced 
to a convergence. We will see this explicitly in Section 4. 

One may ask where the complex structure and Kahler moduli enter in this setup. The 
complex structure enters implicitly through the basis of holomorphic sections s a . As for the 
Kahler class, this is ci(C). Of course, the Ricci flatness condition is scale invariant, so the 
overall scale is irrelevant; however the point of this is that if h l>1 > 1, then by appropriately 
choosing £ we choose a particular ray in the Kahler cone. 

Therefore, if we can find the unique balanced metric for a given C, we have a candidate 
approximation scheme. Let us evaluate it by our criteria. One great advantage is that we 
have a control parameter k, which is easy to implement and is mathematically natural. The 
balanced metric is natural in many respects, which should make it possible to get error 
bounds like Eq. (4), though this has not yet been done. 

On the other hand, the balanced metric is not in general the most accurate approximation 
within this class of metrics to gRF- As discussed in [8], one can find other series of metrics 
which converge to gRF faster than 1/k 2 . 

As an illustration, once we choose our measure of accuracy e, say Eq. (3), we can propose 
a simple scheme which is guaranteed to be the most accurate possible. It is a two-step pro- 
cedure in which we take the balanced metric as the starting point for a numerical search in 
the space of parameters h for the metric which minimizes e. As we discuss a bit later, stan- 
dard numerical optimization routines will work for this purpose if the starting point is close 
enough to the actual minimum. This would not be the most efficient possible approach; one 
could improve the efficiency by using information about the linearized problem, as discussed 
in [8]. 

For the applications we have in mind, for example being able to detect anomalously 
large or small numbers in observables (perhaps having to do with singularities), control, 
mathematical naturalness and ease of programming are more important than accuracy, and 
thus we stick to the balanced metric in our present work. 

As another simple physical illustration, suppose that by using the techniques of [10] we 
could use these results to get canonically normalized fields and physical Yukawa couplings in 
quasi-realistic compactifications. At this point we would probably be much happier to have 
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results for quark and lepton masses in a variety of models which were guaranteed accurate 
to within a factor of 2, than to have "probable" 10% accuracy in one model. Of course this 
is a rather long term goal, but the point should be clear. 



3 Numerical integration on Calabi-Yau varieties 
3.1 Basic setup 

It is clear from the outset that analytic evaluation of the integrals appearing in the T-map 
(7) is not possible. On the other hand, if the integrands are smooth and relatively slowly 
varying functions, it will be possible to evaluate the integrals using Monte Carlo methods. 
This is clear for the sections themselves. Since h is positive definite, the denominator in (7) 
is strictly positive, mitigating (though not eliminating) the possibility of numerical blow-ups. 

Let X be a compact Calabi-Yau n-fold, 3 with its corresponding holomorphic n-form 
Q G A n,0 (X). The volume form Q A Q determines a natural measure d/j,Q on X in the sense 
that 

' fdfi n = [ fQAQ. 
x J x 

From now on we will not distinguish between a top form and the associated measure. 

We can use dfiQ to measure volumes. For an open set U C X the indicator or character- 
istic function 1^ is defined by 



The measure of U is its volume 



1 if x E U 
ifx£U. 



A*n(W) = / ludfin = vol(W). 
Jx 



To do a Monte Carlo integration, one would ideally like to produce samples of points 
on X which are uniformly distributed according to the measure dfi^. This means that for 
every sample of points {qi E X}fj! v the expected number of points within each open subset 
Wclis 

ST 1 („ \ AT I* 1 ^ 

Using this, we can estimate integrals as finite sums: 

r i Np 

/ /^w^-^/fe), (9) 

Jx V i= i 



3 Although most of what we present generalizes to varieties other than Calabi-Yau, we restrict attention 
to these spaces due to their importance in string theory. 
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The statistical error of such an approximation is of order 1 / y^/Vp times a quantity propor- 
tional to the mean of the /(gi)'s. [14]. 

In practice, producing samples of points which are distributed according to the measure 
/jLq is not so easy. One way to overcome this problem is by producing points which are 
uniformly distributed according to another auxiliary measure, say dfiA- Let us assume that 
d/j,A is associated to the global top form A. The ratio Q A Q/A is a global function on X, 
which we call the mass function vtla- At a point x it is defined to be the ratio of the two top 
forms evaluated at x : _ 

m ^> = -a® 1 - < 10) 

In general this function is neither constant nor holomorphic. 

While one could use this information to generate a sample distributed according to dfia 
(e.g., by rejection sampling or MCMC), it is simplest to explicitly put the mass function into 
the integrand. Thus, given a sample of points distributed according to dfiA, and the mass 
function, we can estimate Eq. (1) as 

/ fdto = [ fp^d^A* ^ £ /tow*), (ii) 

Jx Jx dfl A L m j fcf 

The presence of the mass function increases the statistical error. On the other hand, the 
generic values of our mass function are order one, and this is a very mild penalty. 

Rather than regarding the Monte Carlo as a way to estimate the original T-map, an 
alternate point of view is to regard the right-hand side of Eq. (11) as defining a new measure 
v and a new T-map, leading to a new //-balanced metric which approximates the desired 
Q A 0-balanced metric. An advantage of this point of view is that in [8] it is shown that 
(under a very mild hypothesis on v) the new T-map is contracting, and the new ^-balanced 
metric is unique. Thus, numerical pathologies will not enter at this stage, provided that we 
use the same sample of points throughout the computation of the balanced metric. This 
is also advantageous for efficiency reasons, so we always do this. One can then repeat the 
computation with different samples to estimate the statistical error. 



3.2 Generating the sample 

We now discuss how to efficiently generate points according to a known simple distribution. 
In this paper we restrict to the case of X a degree d hypersurface in P n+1 . For definiteness 
let X be defined as the zero locus of the degree d homogenous polynomial /. The case of a 
complete intersection is a straightforward extension. Our main interest will be in d = n + 2, 
but we can be more general for the time being. 

First, it is easy to generate random points distributed according to the Fubini-Study 
measure (for any h) in the ambient P n+1 . We simply generate uniformly distributed points 
on S 2n+3 , a standard numerical problem, and then mod out the overall phase. 

Using this distribution, one approach to generating points on X would be to keep only 
the points that lie sufficiently close to X, in other words satisfy the defining equation of X 



7 



with a given precision, and then use a root finding method (say Newton's method) to "flow" 
down to X. In essence this is a rejection- type algorithm. We implemented this strategy, 
but it has some problems. First, it is hard to control the emerging distribution on X (this 
depends on details of the root finding method). Second, it is an order of magnitude slower 
than the second method we are about to describe. 

The approach we use starts by taking a pair of independently chosen random points 
(X,Y) G p n+1 xP n+1 , which define a random line in P n+1 . By Bezout's theorem, a generic 
complex line in P n+1 intersects X in precisely d points, and we take these d points with equal 
weight. Repeating this process M times generates some random distribution of N p = dM 
points. 

One advantage of this approach is that finding all d roots of f(z) = numerically is 
not much harder than finding one root. But the main advantage, as we show in Section 3.3 
using results by Shiffman and Zelditch on zeroes of random sections, is that that the resulting 
points are distributed precisely according to the Fubini-Study measure restricted to X. The 
mass function (10) is then computable quite efficiently 

A possible disadvantage for some applications is that the resulting sample will have 
correlations between the points in each d-fold subset. For our purpose of Monte Carlo 
integration, this is not a problem, as (11) is the expectation value of a function of a single 
random variable, and does not see these correlations. If one were considering functions of 
several independent random variables, one would probably want to further randomize the 
sequence (say by permuting points between subsets) to remove these correlations. 

3.3 Expected values of currents 

Let us start with a smooth compact algebraic variety X, and an ample line bundle £ on X. 
As reviewed in Section 2, this means that C k defines an embedding i& into projective space 
for any k > ko, for some positive integer fc : 



The idea is to consider random global sections of C k , distributed uniformly according to a 
natural measure, and look at the expected value of the zero locus that they cut out in X. 
For this it is convenient to use the Poincare dual formulation, where the divisor associated 
to a section becomes a form, and ask what is the expected value of the random forms. 
Shiffman and Zelditch answer this question, and the more general one when we intersect I 
such divisors, in full generality using the language of currents. This section is a brief review 
of some aspects of their work [15, 16]. For brevity we adapt their results to fit our needs, 
rather than reproducing them verbatim. 

The space of global sections T = H°(X, C k ) is a complex vector space of dimension d^. If 
we choose a basis for it, then it automatically defines a hermitian inner product, with respect 
to which the basis in question is orthonormal. Conversely, given a hermitian inner product 
(-, •) on T, there is an orthonormal basis B = {s\, . . . , Sd k } on T. Now given s G T, we can 
expand it in the basis £>, and the inner product induces a complex Gaussian probability 




(12) 
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measure on T: 

rf 7(s) = ^e- ||c|l V*c, where s = J2f =1 c jSj and ||c|| 2 = Z%i \cj\ 2 " ( 13 ) 

Given a metric h on the line bundle C, as explained in (6), h defines a hermitian metric on 
T. This is the inner product that we are going to use on T = H°(X, C k ) throughout this 
section. 

Given a random variable Y on the probability space (r, dj), the expected value of Y in 
the probability measure is 

E(Y) = J Ydj. (14) 

We can think of the probability space (r, d'-f) in a slightly different way. Consider the 
unit sphere 

ST = SR°(X,C k ) = {s e R°(X,C k )- (s,s) = 1}. 
The Gaussian probability measure on T restricts to the uniform measure d[i on ST. The 
expected value of Y\g r is E(Y) = J ST Yd\i. On the other hand, the uniform measure on the 
sphere ST descends to the Fubini-Study measure on the projectivization FT. This alternative 
view will be very useful later on. 

If we choose a section s E T = H°(X, C k ), then there is a divisor Z s associated to it, 
which, roughly speaking, is the zeros of s minus the poles of s. Since we work with L k very 
ample, Z s consists of only the zero locus of s. Given the probability measure d^ on T, we can 
choose s randomly, and ask what is the expected value of the random variable Z (defined 
by s t— > Z s ). This same question can be asked in an equivalent form using Poincare duality. 
The Poincare dual of Z s is a (1, 1) form T s , and it is more convenient to work with forms in 
this context than to work with divisors. In general T s is not a C°°-form on X, but it can be 
given an explicit expression using the notion of currents, i.e., distribution valued forms. 

Currents are defined as it is customary in the theory of distributions. 4 Let Qq' 9 (X) 
be the space of compactly supported C°° (p, g)-forms on X, and for now we assume that 
dimX = n. The space of (p, q)-currents is the distributional dual of Vt^~ p ' n ~ q {X): T> p,q (X) = 
Q™-p> n -i(Xy . An element of T> p,q (X) is a linear functional on VL^~ p ' n ~ q (X) which continuous 
in the C°° norm. 

The usefulness of currents in our context stems from the fact that Poincare dual Ty of 
an algebraic subvariety Y, defined by 

/ T Y A a = / L*(a), for any a G QQ imY,dimY (X) , 
Jx Jy 

oftentimes has an explicit form in terms of currents (l: Y <^-> X is the embedding). Let us 
focus on the case when Y is a hypersurface, or more generally the zero divisor of a section 
s G T = B°(X,C k ). In this case the current is given by the Poincare- Lelong formula: 

T s = -dd log(s,s) G V l ^(X). 

IX 

For a nice introduction to distributions and currents in algebraic geometry the reader can consult [17], 
Chapter resp. Chapter 3. 
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T s is also known as the zero current of s. Thus, the Poincare-Lelong formula induces a map 

T:PH°(X, C k ) — >V 1 ' 1 (X), s^T s . 

As discussed earlier, the Fubini-Study measure makes PH°(X, C k ) into a probability 
space, and we can also view T as a random variable. Since the currents form a linear space, 
we can inquire about their expected value in this probability measure 

E{T)= [ T s dfx FS (s), 

JPH°(X,C k ) 

As it happens oftentimes in the theory of distributions, although T s is not a C°° form, E(T) 
is, and we have the following proposition: 

Proposition 3.1 ([15, Lemma 3.1]) With X and £ as above, for k sufficiently large so that 
Kodaira's map associated to H°(X, C k ), as defined in (12), is an embedding, the expected 
value of the random variable T representing the zero current is 

E{T s ) = ±ilu F k s , 

where u k s is the Fubini-Study 2-form on PH°(X, C k ), and i* k is pullback of forms (in this 
case restriction). 

This result generalizes to the case when we intersect several divisors, and this will be 
the case of main interest to us. Since intersection of subvarieties is Poincare dual to the 
wedge product, there is an obvious guess how Prop. 3.1 should generalize. Let s±, . . . ,s m be 
sections of H°(X, C k ), and consider the zero current Tj™^ Sm associated to the set 

Z S i,...,s m = {x G X: si(z) = • • • = s m (x) = 0}. 

It is obvious that if we consider any m linear combinations of si, . . . , s m (which are themselves 
linearly independent), then they determine the same zero set Z Su ___ jSm . Hence Z sli ... |Sm is really 
associated to the m-plane spanned by s±, . . . , s m in H°(X, C k ). Therefore the probability 
space is the Grassmannian of m-dimensional subspaces of H°(X, C k ), with its natural Haar 
measure d[iu aar (a generalization of the Fubini-Study measure). So we can ask what is the 
expected value of this zero current, computed with dfiHaar- 

Proposition 3.2 ([15, Lemma 4-3]) In the notation of Prop. 3.1, the expected value of 
the zero current T^ m \ associated to the simultaneous vanishing of m random sections of 
H°(X, C k ), distributed according the Haar measure of the corresponding Grassmannian, is 

E (T {rn) ) = k™- 1 (i* k ul s ) m (15) 
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Note in particular, that unlike in Prop. 3.1, the distribution is according to the Haar 
measure, but the final result still involves the Fubini-Study form on H (X, C k ). This is 
in fact natural, given that an m-tuple of sections, each distributed according to Fubini- 
Study on E°(X, C k ), is the same thing as one m-plane, distributed according to Haar on the 
Grassmannian. 

As usual, besides having expected values, random variables also have variance. The zero 
current is no different in this respect. Its variance has been computed recently in [16, 
Theorem 1.1]. In particular, it was shown that the ratio of the variance and the expected 
value goes to zero as k is increased 



[Var (T( m ))] 



1/2 

m 1 

~ k, 2 4 



s(rw) 

3.4 The Calabi-Yau case 

In this section we put all the pieces together, and explicitly show how to build the numerical 
measure {qi G X, m(qi)} i ^ 1 introduced in Section 3.1. We will focus on a smooth Calabi-Yau 
hypersurface X in P n+1 of degree n + 2. Let X be given by the zero locus of the degree 
n + 2 homogeneous polynomial /, and let (Zq, Z\, . . . Z n+ i) be the homogeneous coordinates 
on P n+1 . We denote the embedding by 

i: X = Z(f) ^ P n+1 . 

Our approach is to generate random points on X using random lines on P n+1 , by look- 
ing at the intersection of these random lines with X. We can view a random line as the 
intersection of n random hyperplanes. This allows us to compute the expected value of the 
corresponding zero current using the techniques of Section 3.3. 

For computational purposes, designing an algorithm to generate points on X in such 
a fashion is straightforward. To generate a random line on P n+1 we generate two random 
points, which lie on the unit sphere S 2n+3 C C n+2 , and are distributed uniformly on this 
sphere. For instance, to generate random points uniformly on S 2n+3 we can start with 
the unit cube in M 2n+4 , i.e., [—1, l] 2n+4 c M 2n+4 . Using a good quality random number 
generator we generate an uniform distribution of points in [—1, l] 2n+4 . Now take only those 
points which fall within the unit disk D 2n+4 , and then project them radially to the boundary 

The intersection of the random line with X can be computed by restricting the defining 
polynomial / to the line. As a result, computing the common zero locus reduces to solving for 
the roots of a polynomial of degree n + 2 in one variable. We find numerically the n + 2 roots 
using the Durand-Kerner algorithm [18, 19], which is a refinement of the multidimensional 
Newton's method applied to a polynomial. This whole approach turns out to be very efficient 
in practice, in that one can generate a million points on a quintic in a matter of seconds. 
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3.4.1 The expected zero current 

We chose to work with the hyperplane line bundle C = OxiX) on X. OxiX) is ample, and 
its global sections are in one to one correspondence with the homogeneous coordinates of the 
ambient P n+1 (restricted to X). The associated Kodaira embedding is precisely the defining 
one: 

h:X^ P n+1 = PH°(X, O x {\))\ 

If we take n sections of C = OxX)i and look at their common zero locus, then by Bezout's 
theorem this is generically n + 2 points (degenerations might occur). Therefore, considering 
random n-tuples of sections will give random n + 2-tuples of points on X. But now we can 
tell how these points are distributed, provided that the sections were distributed according 
to the Fubini-Study measure on PH°(X, OxiX)) = pn+1 - Usin S Pr °P- 3 -2 we know that the 
expected value of the zero current associated to the n + 2 points of intersection is (ij u>^£+i) . 
This is an (n,n) form on X, and plugging it into (10) we obtain the mass formula 

m ( x ) = (■* fs v* ( z )- ( 16 ) 
\i 1 u; pn+ i ) 

3.4.2 The numerical mass 

Let us look at the two differential forms involved in Eq. (16). For this we first choose affine 
coordinates w a = Z a /Z , % — 1, 2, . . . , n + 1 on P n+1 . The Fubini-Study 2-form on P n+1 is 



U! 



pn + 1 

FS 



dw a A dw a ( Yl w adw a ) A ( u! a dw a ) 



(17) 



The pullback (il uip£ +1 ) n £ f2g' n (X, Z) is a top form on X. Let X\, . . . , x n be local coordinates 
on X. Then 

(a* , ,FS \ _ ® W a (, ,FS \ d W b f-,o\ 

and (ij is proportional to the determinant of this matrix. For obvious reasons we 

need to evaluate this determinant. Let us outline how this can be done, paying attention to 
some of the numerical aspects. 

The idea is to choose local coordinates on X that are convenient to work with. Let us 
start with the point P on X with homogenous coordinates Zi. To minimize the numerical 
error we go to the affine patch where \Zi\ is maximal. Without loss of generality let us 
assume that this happens for i = 0. The affine coordinates are w a = Z a /Z . 

Let p be the affine form of/, i.e., p(w) = /(l, w 2 , ■ ■ ■ , w n+ i). This equation determines 
one of the w a -s in term of the others, as an implicit function. Let us assume for the sake 
of this presentation that dp/dw n+ \(P) ^ 0. The implicit function theorem then tell us 
that in an open neighborhood of P w n+ \ is a function of the remaining variables: w n+ \ = 
w n+ i(wi,W2, ■ ■ ■ ,w n ). This allows us to choose the coordinates u>i, . . . ,w n to be the local 
coordinates x±, . . . , x n on X. 
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This choice of coordinates is quite advantageous for computing (18). All we need is to 
compute dw n+ i/dxi, as dwj/dxi = 5y. This can be done algebraically, without explicitly 
solving the p = equation. Namely, using the fact that 

p(w ly ...,w n , w n+1 (w!, ...,w n )) = 

is the identically zero function, its derivative with respect to any Wi vanishes identically, for 
i — 1, . . . , n. As a result we have that 

For numerical stability one should always solve for the variable for which \dp/dwi(P) \ is the 
largest. 

The second differential form entering (16) is Q A Q. The holomorphic n-form Q £ 
Q n '°(X, C) can be represented using the Poincare residue map [17, Section 1.1] 

n _ ( i y-i du?1 A dw2 • • • A dwj A . . . A dw n 
~ [ ' dp(w)/dwi ' 1 } 

where dwi means the omission of dwi in the wedge product. 

These explicit expressions allow us to perform integrals numerically on elliptic curves, 
K3 surfaces and more interestingly, quintic 3-folds. 



3.4.3 Symmetries 

Suppose our Calabi-Yau X is preserved as a complex manifold by the action of a discrete 
group T. Then a Ricci flat metric whose Kahler class uo is preserved by V will also be T- 
invariant, because it is unique. As we will see in this section, the same statement applies to 
the balanced metrics as well. 

A general hermitian N by N matrix has N 2 independent real coefficients. On the other 
hand, if the Calabi-Yau X has discrete symmetries, then we expect to find symmetry relations 
between the matrix elements of T(h). Taking advantage of these relations can drastically 
reduce the size of the problem. In this section, we argue that these symmetry relations are 
respected by the balanced metric and the T-map, and explain how we used them in the 
examples of Section 4. 

Next, let us review the symmetries of X defined as a hypersurface in P n+1 by the degree 
n + 2 homogenous polynomial 

n+1 n+1 

/ = E z r +2 -^+ 2 )^n^ ( 21 ) 

i=0 i=0 

Here ip controls the complex structure of the hypersurface. Using the fact that X is Calabi- 
Yau, the symmetry group is finite. To find the symmetries of X we consider two natural 
group actions on P n+1 , and impose conditions such that these group actions descend to X. 
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We start with the abelian group 

n+1 

0Z P C GL{n + 2) 

i=0 

that acts by independently rescaling the n + 2 homogenous coordinates i— > ajZj, where 
aj are pth roots of unity. Since the projective coordinates are defined only up to overall 
rescaling, we have to mod out by the diagonal action AZ p and find the group 

n+1 

Z P /AZ P C FGL(n + 2) 

8=0 

acting on P n+1 . In order for this group to descend to X, it must leave the defining equation 
(21) invariant. In the Fermat case, that is ip = 0, we set p = n + 2 and find that the 
Calabi-Yau is invariant under 

n+1 

Z n+2 /AZ n+2 = (Z„ +2 )" +1 . (22) 

i=0 

For non- vanishing ift, the «j have to obey the additional constraint n^o 1 tt « = 1- This shows 
that the symmetry group is a subgroup of Z n+2 /AZ n+2 — (Z n+2 ) n+1 , given by the kernel of 
the product map (ao, . . . , a n +i) ^ YYi=o a i- We call this group Ab n+ 2, and it is clear that 
there is an isomorphism Ab n+ 2 = (Z n+2 ) n . For example, in the case of the torus defined by 
our cubic in P 2 , Ab% = Z3. At the Fermat point this group is enhanced to Z|. 

The second symmetry group we consider is the symmetric group on n + 2 elements § n +2 ■ 
This group acts by permuting the coordinates of P n+1 . Since Eq. (21) is invariant under 
permutations, S n+2 is a symmetry of X as well. 

To see how these actions on the coordinates of P n+1 induce an action on {s a }, the global 
sections of the line bundle L k on X defining the embedding in P^ -1 , we can use some simple 
algebraic geometry. The fact that X is given by a hyperplane in P n+1 gives a natural way 
to parameterize the global sections of C k = Ox(k). We start with the short exact sequence 
(SES) defining X: 

^(9pn+i(— n — 2)— ^C P n + i ^O x ^0 

Tensoring with 0p»+i(fc), and using the fact that H^P™ 4 " 1 , ¥ n+i(k - n - 2)) = 0, for k > 0, 
we get another SES: 

H°(P n+1 , ¥ n+i(k -n- 2))—^ H°(P n+1 , Opn+i (jfe)) H°(X, C k ) -0 (23) 

which shows that the global sections of H°(X, C k ) can be parameterized by degree k mono- 
mials in n + 2 variables modulo the ideal generated by /. Therefore the sections inherit an 
obvious group action. 
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In particular, we also find that 



N = dimH°(X,£ fe ) = 
In addition, note that the map : X ^ P^ 1 factorizes 



n + k + l\ f k-1 
k J ~ \k-n-2j' 




(24) 



The second embedding, v: P n+1 <^-> P 7 ^ 1 , is the (Veronese) map associated to the incomplete 
linear system on P n+1 induced by the complete linear system \C k \ on X. 

We will now consider the consequences of these actions on the T-map and the sequence 
of hermitian matrices {T'(/i)}; =12i .... First we consider the action of Ab n+ 2- We assume 
that T°(h) is invariant under the group action. This is a choice we can always make. Since 
Ab n+2 is an abelian group, its irreducible representations are one dimensional, and can be 
labeled by the characters. Each section s a transforms under a character Xa- The operator 
T is defined in terms of the sections, and the Ab n+2 will force some of the T(h) a p matrix 
elements to be zero. To better understand this we look at a toy example: the integral of an 
odd function a on R. The group G in question is Z2, and acts on R by x 1— > —x. Z2 has only 
one nontrivial representation, and being odd, a transform in this irrep. Now we have that 



a(t) dt = a(—x) d{—x) — — I a(x) dx 



-DC 



where we used a change of variable t = —x. This implies that f*™ a{x)dx = 0. 

More generally, in R n for a function a, a group G, and an element g £ G, we can do the 
change of coordinates t = g ■ x and then 



a(t)dV(t)= / a(g-x)g*dV(x) = / a(g • x) dV{x) . (25) 

X JgX JX 

Here we assumed the measure to be G- invariant. 

Applying Eq. (25) for G = Ab n+ 2, and using the fact that s a transform as a character of 
G = Ab n+2 , it gives that 



T{h) a - p = uv , rrrTz dfi n = Xa{u)Xp(u) T{h) a - p . (26) 

VOl^vY ) J x \\ s \\h 

We used the fact that X is invariant under that action of Ab n+2 i and so is the measure Q Afl, 
and the denominator (The latter follows by induction from the initial choice of T°(h) 

being invariant.) In particular, if Xa(u)xp(u) 7^ 1, for any u £ Ab n+2 , then the corresponding 
T(h) a p has to vanish. In our numerical routine we impose this vanishing condition on all 
the matrices T l (h). 
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A similar argument applies for G = § n +2- Since §„+2 is not abelian, and hence its 
generic irreducible representations are not one dimensional, this constraint does not result 
in vanishing rules, but rather sets a priori independent coefficients of T(h) equal to each 
other. To see how § n +2 acts, recall from Eq. (23) that 

H°(X,£) = H°(P n+1 ,£Wi(l)) = C n+2 

is the fundamental representation of S n +2, call it F. Then H (P n+1 , Ofn+i (k)) is the /cth 
symmetric tensor power of F, Sym k F, and by (23) H (X, C k ) is a quotient of this. Now we 
can return to Eq. (25). Once again, we choose T°(h) to be invariant under § n +2, and then 
induction and Eq. (25) tell us which matrix elements of T(h) equal each other. 

Therefore, imposing the symmetries of both finite groups, the number of independent 
components of T l (h) (for any /) reduces significantly. To illustrate this we consider k — 12 
on the quintic in P 4 , i.e., n = 3 (this was the largest k we computed). In this case N = 1490. 
This means that T(h) is a hermitian matrix with 2,220,100 components. Taking into account 
the Ab§ and §5 relations, one is reduced to computing 9800 components. This simplification 
speaks for itself. 



4 Numerical results 

In this section we present our explicit numerical results. The main object that we compute is 
the balanced metric associated to the embedding of the quintic threefold defined in Eq. (21). 
For definiteness we chose to work with ip = 0.1, but also tested other values of ip. We also 
considered the case of elliptic curves (n = 1) and K3 surfaces (n = 2). In all these cases we 
obtained results similar to the ones to be presented here. 

To find the balanced metric we study the associated N x N matrix hk for several values 
of k, from k — 1 to k — 12. We use hk to construct the associated Kahler form Uk on X, 
and check how well it approximates the Ricci flat metric. We do this in several ways. 

First, one can study the function defined in Eq. (2) 

det u)h „ _ 
ri k = S : X — > R. 

For a good approximation g k to the Ricci flat metric gup the function is almost constant. 
We study the behavior of statistically, by summing over all the regions of X, and also 
locally paying attention to certain special regions of the threefold. 

Second, one can compute the Ricci tensor of Uk- To check pointwise how close to zero 
the Ricci tensor is, we need a diffeomorphism invariant quantity. We chose to work with the 
Ricci scalar. We also perform this analysis for several values of k, and show how the Ricci 
scalars decrease pointwise with k. 

Before presenting the results let us comment on the errors coming from Monte Carlo 
integration. We estimate them by computing the balanced metrics associated to different 
samples of points, and then looking at the mean and variance of each individual matrix 



16 



element. Ideally, one would like to produce samples of points with minimal induced error. 
Constructions that reduce the standard deviation of the integrals are refinements to the 
theory of numerical integration presented here. Markov Chain Monte Carlo techniques, con- 
struction of lattices on Calabi-Yau varieties, and of quasi-random points on such manifolds 
are different approaches that one could consider. 



4.1 Approximating volumes v.s. Calabi-Yau volume 

Here, we consider the way the function 

\r} k -l x \'X — >R+, x (->■ \r}k{x) - 1| 

behaves on X. As argued earlier, we expect \r)k — lx I to approach the constant zero function. 
One can study the deviation of \r\^ — \x\ from the zero function by computing the integral 



Jx 



(27) 



We compute this integral by our Monte Carlo method, which introduces an error, and this 
error can be estimated by 



5 a i 



(\Vk - lx| - erf dfi S 



1/2 



where N p is the number of points used to perform the Monte Carlo integration in (27). 



(28) 
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Figure 1: au and Ricci scalars. 

In Fig. 1 we plot the values defined in (27) for k — 3, . . . , 12. The error bars for each 
value are the corresponding standard deviations (28). We also see how the errors decrease, 
along with <Tfc, for higher and higher k. The fit in Fig. 1 is a curve of type 



17 



as we expect from the theory. 

We can also study the local behavior of r]k by restricting it to a subspace. Given our 
quintic 3-fold, we consider the rational curve defined by 

(Zq = zq, Z\ = —zq, Zi = zi, Z3 = 0, Z4 = —zi), (29) 

where Z t are homogeneous coordinates on P 4 , while (z , zi) are homogeneous coordinates 
for P 1 . This rational curve lies on every quintic defined by Eq. (21). 




Figure 2: The values of r\ on the rational curve, for k — 1, 3, 4, 5, 7, 9, 11 and 12. 

In Fig. 2 we plot the values of function 77^ restricted to the rational curve defined 
above for 12 different values of k, ranging between 1 and 12. More concretely, given the 
embedding (29), we choose the local coordinate system on P 1 defined by t = z\j 'zq, and take 
the stereographic projection of the t-plane. Using spherical coordinates (9, <ft) on P 1 = S 2 
we embed it into P 3 , by the parameterization 

Zq = sin 9 cos <p, Zi = sin 9 sin + i cos 9. 

In the radial direction of M 3 we plot the function 77^. As expected, 77^ approaches the constant 
function 1 as k increases. 

4.2 Ricci scalars 

Next we discuss how we compute the Ricci scalar on X. The Kahler potential on X is given 
by the restriction of the Kahler potential (5), which is associated to a balanced metric. The 
metric and Ricci curvature on X are given by 

g h = ddK h , Ric = dd det g h . (30) 
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To compute these quantities we need the first and second derivatives of the sections s a 
with respect to the coordinates of X. We can compute these derivatives algebraically, using 
the ideas outlined in Section 3.4.2. To get the actual value of g and Ric, at a specific point 
x G X, we evaluate our numerical expressions at this point. Let us stress that this way g 
and Ric are evaluated algebraically, without numerical derivatives. 

Let us now discuss the effect of the parameter k on the Ricci scalar R. The theoretical 
estimate predicts a vanishing of R for large k as 

\\R k \\<l + ^ + ---, (31) 

where 7 and 5 are constants, in any C r norm. We can see this pointwise, after a short 
statistical analysis. Let us pick 100 randomly chosen points Pi on the 3-fold, and compute 
the associated sets of Ricci scalars R k , for k G [3, . . . , 10]. To obtain a normalized set R k of 
Ricci scalars we rescale all R k by I/R3, and we do this for every point. This normalization 
leads to a more meaningful comparison between the different points. 

We are only interested in points where R k shows a generic behaviour. For this we compute 
the expectation value of |i?io| {k = 10 gives a good accuracy) 

1 100 

<l^°l> = Too^ lRl ° m = - 48 - 
1=1 

We consider a point P generic if it lies at a distance of order one from the mean (|i?io|), that 
is 

|£ 10 (P)-(|i2 10 |)|<0.5. 

We find that 95 points out of 100 obey this condition. We use these 95 points for our 
statistical check of Eq. (31). In Fig. 1 we plot (|-Rfe|) as a function of k. We do a least square 
fit for both 7 and S, and obtain very good agreement. 




Figure 3: Histogram of Ricci scalars, and scattered plot of mass versus the Ricci scalars. 
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We can obtain a visual picture of how the Ricci scalars decrease with k by plotting them 
on a histogram. We use the same 100 points from above. To keep the picture simple we plot 
only six k values: 3, 4, 5, 7, 9 and 12. Unlike earlier, here we do not normalize the Ricci 
scalars pointwise, instead we reorder them decreasingly. This ordering is done for no other 
reason but to enhance the visual clarity. Accordingly, in the first graph of Fig. 3 we plotted 
the value of the Ricci scalar for every point (the y axis is the absolute value of Ricci scalar, 
while on the x axis we have the points from 1 to 100). We did this for the k values indicated 
above, and used different colors to distinguish them (e.g., yellow corresponds to k — 3, while 
k = 12 is blue). It is evident from this graph that by going from k = 3 to k = 12 the Ricci 
scalars decrease by an order of magnitude, in line with the theoretical expectation. 

One observes that for any k there are a few anomalously large Ricci scalars. To un- 
derstand this we can do a scattered plot of the mass of that point versus the Ricci scalar 
in question (we present this for k = 12), using 1000 random points. This is the second 
graph in Fig. 3. The picture shows a correlation between large values of the Ricci scalar and 
large values of the mass (large in a logarithmic sense). In other words, in regions where our 
point generator needs large correction, via the mass, the balanced metric is a less accurate 
approximation of the Ricci flat metric, compared to points with smaller mass. This fact is 
then amplified by the formula for the Ricci tensor in Eq. (30), where the logarithmic scale 
is also supplied. 

Finally let us note that l/r/i is precisely the mass function (10). This is because the k = 1 
balanced matrix is proportional to the identity, a consequence of the discrete symmetries 
present for our quintic. This leads to an alternative interpretation for the first graph (k = 1) 
in Fig. 2, as depicting the inverse masses of the points on that rational curve. 

4.3 Discussion 

We will discuss further applications of these results elsewhere; here we discuss the advantages 
and limitations of this approach compared to others, for example position space methods 
[20]. 

The runtime of a computation of the balanced metric can be approximated as 

T = N it x N p x S 2 , 

where S is the number of independent sections (taking into account discrete symmetry), 
N p is the number of points used in the Monte Carlo integration, and N it is the number 
of iterations of the T-map required for convergence. Since convergence is exponential, this 
leads to a rough scaling with the accuracy as 

T ~ ^S*. 

The value of S required for a given accuracy depends on the symmetries and dimension. 
For the balanced metrics, we expect to need k ~ 1/ \ft\ as discussed in section 2 this 
could probably be improved by choosing a different scheme if accuracy were paramount. 
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For hypersurfaces in n complex dimensions, we then have S ~ N ~ k n+1 , leading to a 
rough overall scaling of T ~ l/e n+3 . This might be compared with a (naive) T ~ l/e 2n 
for position space methods, so the two appear generally competitive. However, along with 
the other advantages we mentioned, we believe the approach we are discussing is far easier 
to program, and requires relatively little effort to adapt to different manifolds, and related 
problems such as hermitian Yang-Mills. 

Since the sections s a of Ox{k) are degree k polynomials, this basis is a simple type of 
Fourier or momentum space basis. Very roughly speaking, a degree k basis should be able 
to represent arbitrary structures on length scales down to 1/k. They are particularly well 
suited for approximating smooth functions, as the Fourier coefficients of such a function 
fall off faster than any power of k (see the appendix of [8] for more precise statements). 
This is advantageous as the Ricci flat metric is smooth, suggesting that other approximation 
schemes could do better than e ~ 1/k 2 . 

On the other hand, in some limits (say a conifold limit) the metric can develop structure 
on small scales, which might not be well represented by a fixed k basis. This is also a problem 
for position space methods with a fixed lattice; there one deals with it by multi-scale methods, 
for example allowing the lattice spacing and structure to vary over the manifold. This is very 
powerful but also very intricate to program. In the present context, rather than increase k, 
one might look for analogous simplifications; either a multi-scale method which uses different 
k in different regions (or even some sort of wavelet-inspired method). Or, since we have many 
explicit expressions for Ricci flat metrics near singularities, it might be useful to develop a 
way to patch these solutions into the global approximate solutions we discussed. 

On the computer code 

Our numerics is based on code that has been written entirely in C++. Our experience shows 
that these computations must be done in a compiled language, rather than an interpreted 
one. We have made extensive use of the following Boost libraries: uBlas, random, bind and 
thread. These libraries are on par with Fortran code, due to implementation techniques 
using expression templates and template metaprograms. The computations were done on 
an Athlon 64 4800+ dual core machine, with 4GB memory. The computational time ranges 
from minutes, for low k, to hours, and eventually 2 days (for k = 12). 
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